Three-Body Bound States in a Lattice 
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We pursue three-body bound states in a one-dimensional tight-binding lattice described by the 
Bose-Hubbard model with strong on-site interaction. Apart from the simple strongly-bound "trimer" 
state corresponding to all three particles occupying the same lattice site, we find two novel kinds of 
weakly-bound trimers with energies below and above the continuum of scattering states of a single 
particle ("monomer") and a bound particle pair ("dimer"). The corresponding binding mechanism 
can be inferred from an effective Hamiltonian in the strong-coupling regime which contains an 
exchange interaction between the monomer and dimer. In the limit of very strong on-site interaction, 
the exchange-bound trimers attain a universal value of the binding energy. These phenomena can 
be observed with cold atoms in optical lattices. 
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Introduction. For decades, many-body physics in pe- 
riodic potentials has been one of the central topics of 
research in condensed-matter theory The complex- 
ity of such systems requires the use of model Hamil- 
tonians which, for sufficiently deep lattices, take into 
account only short-range interactions and neglect cou- 
pling between different bands of the periodic potential 
2j. Despite apparent simplicity, tight- binding many- 
body Hamiltonians, such as the Hubbard model, have 
been able to elucidate and predict many important lat- 
tice phenomena. A prime example is the transition from 
a superfiuid to a Mott insulator phase [H experimentally 
realized with cold bosonic atoms in optical lattices Q • In 
another recent experiment [f|, stable repulsively-bound 
pairs of atoms in an optical lattice were observed. This 
seminal achievement has led to several theoretical stud- 
ies of the dynamics of atom pairs in periodic potentials 
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Two decades ago, Mattis [14| established an equation 
for the bound states of three bosons, and predicted Efi- 
mov states [l5[ in a 3D lattice. Here we study the 
three-body problem in a one-dimensional (ID) tight- 
binding lattice described by the Bose-Hubbard model 
with nearest-neighbour hopping and strong on-site in- 
teraction. We find that, in general, for each value of to- 
tal quasi-momentum there are three distinct three-body 
bound states. The first strongly-bound "trimer" state is 
simple, corresponding to all three particles being at the 
same lattice site. The other two hitherto unknown bound 
states arise due to an effective exchange interaction be- 
tween a bound pair — "dimer" — and a single particle — 
"monomer" . These weakly-bound trimers have energies 
slightly below and above the two-body continuum of the 
monomer and dimer scattering states. All three trimer 
states exist irrespectively of the sign of on-site interac- 
tion, while for very strong interaction, the exchange- 
bound trimers attain a universal value of the binding 
energy equal to half the single particle hopping rate. 



Mathematical formalism. We consider bosonic parti- 
cles in a ID lattice described by the Hubbard Hamilto- 
nian 



H = -J 
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where St (bj) is the particles creation (annihilation) op- 
erator and hj = bjbj the number operator at jth lattice 
site, J(> 0) is the hopping rate, and U is the on-site in- 
teraction (which can be attractive or repulsive). We seek 
three particle bound states in momentum representation, 
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where the wavefunction ip(k±, k 2l k 3 ) is symmetric with 
respect to exchange of any pair of particles, with each 
particle quasi-momentum kj G f2 restricted to the 
first Brillouin zone f2 = [— 7r, it]. From the stationary 
Schrodinger equation H = E using the conserva- 
tion of total quasi-momentum K — ki + k 2 -\-k 3 (mod 2tt) , 
we obtain 



tp(ki,k 2 ,k 3 



Mjh) + M{k 2 ) + M(k 3 ) 
' e(fci) + e(fe) + e(*3) - E' 
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where e(k) = — 2Jcos(fc) is the single-particle energy 
(Bloch) band of Hamiltonian (fTJ) , while the functions 
M(k) satisfy the ID Mattis integral equation [14j | 
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M(q) 



e(k) + e(q) + e(K - k - q) - E' 
with 7e(A;) being a generalized Watson integral [l6| 
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I E (k) = - j dq 



e(k) + e(q) + e(K - k - q) - E 
sgn[£ - e(k)]U 
y/[E - e(fc)] 2 - 16 J 2 cos 2 [(if - k)/2] ' 
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FIG. 1: (a) Full three-particle energy spectrum of Hamilto- 
nian ([l]) with U = —10 J, versus the total quasi- momentum 
K. All bound states are obtained via exact numerical solution 
of Eq. Q. (a') Magnified part of the spectrum corresponding 
to dimer-monomer states, (b) Binding energies Eb for the 
off-site (weakly-bound) trimers at K — versus the interac- 
tion strength U < 0. (c) Dimer-monomer spectrum of the 
effective Hamiltonian ([5} • The two bound states are obtained 
via numerical solution of Eq. (|6]). 



Equation (j4|) can be cast as a homogeneous Fredholm 
equation of the second kind with eigenvalue A = 1. 
Hence, for a given U/J and fixed K , it is a nonlinear 
equation for energy E, which can be solved numerically. 
Note also that Eq. ((4| does not imply M(k) to be an 
even or odd function; in fact, there is no symmetry with 
respect to k = 0, unless \K\ = or ir. 

Three-body spectrum. As is well-known 14j . the Bose- 



Hubbard Hamiltonian JTJ) with \U\/J^$> 1 has a narrow 
band of on-site bound states, corresponding to A" > 1 
particles tightly co-localized on the same lattice site, 
with the energy Eb ~ \U(N 2 — N). In the present 
case of A" = 3, we have Eb 3(7. The other two 
known parts of the spectrum of Eq. ([TJ are: the three- 
body scattering continuum of three (asymptotically) free 
particles, with the energy given by the sum of single- 
particle bands, E c3 = e(fci) + e(k 2 ) + e(K — ki — k 2 ); 
and the two-body scattering continuum of a bound pair 
(dimer) and a free particle (monomer), with energy E c2 = 
sgn(U)^U 2 + [4Jcos(Q/2)] 2 -2Jcos(A:-Q), where the 
first term is the energy of a dimer with quasi-momentum 

q BSEl. 

The full spectrum of Hamiltonian (JTJ is shown in 
Fig. QJa). For concreteness, here we consider attractive 
interaction, U < 0, but note that our results equally ap- 
ply to the case of repulsive interaction, U > [17|. As 
seen in the figure, lattice bosons can form two new kinds 
of three-body bound states whose energies E ai and E a2 
lie below and above the two-body continuum E c2 . Some 
properties of these states can be deduced by energy con- 
siderations. First, these are not on-site bound states, 
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FIG. 2: Top panel: Quasi-momentum distributions 
\t/)(kx, fo, K — ki — k2)\ 2 of the three-particle bound states for 
U = —10 J and K — 0; (al) and (a2) correspond to the off-site 
(weakly-bound) trimers with energies below and above the 
two body dimer-monomer continuum, while (b) corresponds 
to the on-site (strongly-bound) trimer. Lower panel: Reduced 
single-particle momentum distributions po(k) for the cases of 
(al), (a2) and (b). 



since their energies E a u2) — U + O(J) are far from 
3U. Next, their binding energies, with respect to the 
E c2 band, are E m ^ 2 ) ~ T^/2 [see Fig. Ufa'), (b)], which 
suggests that these are off-site weakly-bound states of 
a dimer and a monomer. Note that the state above the 
two-body continuum is bound stronger than the state be- 
low the continuum. Finally, they are not Efimov states 
which can exist only in 3D systems near two-boson reso- 
nance [HI Ell EH > since we examine a ID lattice with no 
two-particle scattering resonances (lTj| . 

In Fig. [2] we show the quasi-momentum distributions 
\ip(k\, k 2 , K— k% — k 2 )\ 2 for the three bound states at total 
quasi-momentum K = 0, as well as the corresponding 
reduced single-particle momentum distributions px (k) = 
f n dq\il>(k,q, K — k — q)\ 2 . In cold atom experiments [20(, 
PK{k) can be directly mapped out by absorption imaging 
of the atomic density distribution after sudden release 
from the lattice potential followed by a time-of- flight. For 
the trimer state with energy E a \ lying below the two- 
body continuum E c2 , the Mattis function M(k) at K = 
is symmetric, and the corresponding quasi-momentum 
distribution is peaked at k = 0, similar to that of an 
attractively bound pair of particles @, 0] • For the bound 
state with energy E a2 above the continuum, the K = 
Mattis function is antisymmetric with respect to k = 
0; the quasi-momentum distribution is peaked at \k\ < 
7r resembling that of a repulsively bound particle pair 
0, S 0] • Finally, the on-site bound state (ground state 
for U < 0) has quite a flat momentum distribution, since 
the very large binding energy Eb — E c2 ~ 2U makes the 
on-site trimer practically immobile. 

As can be seen from Fig. [ljb), where we plot the 
binding energies Eb\< 2 ) at K = 0, there are thresh- 
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olds for the existence of full bands of the off-site bound 
states. For the trimer below the two-body continuum, 
the binding energy vanishes when \U\ ~ 4 J: at this 
critical value of U the trimer energy E a i approaches 
the edge of the dimer-monomer scattering continuum 
E c2 = —^/U 2 + 16J 2 — 2 J resulting in K — scattering 
resonance. On the other hand, the trimer above the two- 
body continuum ceases to exist already for \U\ w 8.5J, 
since then its energy E a2 approaches the bottom of the 
three-body continuum E c z — 3e(0) = —6 J (the two con- 
tinua, E C 2 and E c s, overlap for \U\ < 8J). Thus, at 
K = 0, the trimer state with energy E a2 starts to appear 
well in the strong interaction regime, while for larger K 
the threshold is smaller: \U\ ~ 4 J for \K\ — > tt. 

Effective model. The intriguing properties of the off- 
site trimers can be understood in term of an effective 
perturbative model valid in the strong interaction regime 
(\U\/J > 8). Upon adiabatically eliminating the on-site 
trimer states |3j) and states |lj-lfclj) with all three parti- 
cles at different lattice sites, for a single dimer-monomer 
system we obtain an effective Hamiltonian to second or- 
der in J, which has the form 



H, 



H 1 +H 2 + H u 



(5) 



Here Hi — — JJ2j(bjbj+i +H.c.) describes the monomer 
with the single-particle spectrum e(k). Next, 



H, =£< 



is the Hamiltonian for a dimer, with cj being the 

dimer creation (annihilation) operator and rhj — CjCj 
the number operator at site j. The first term in H 2 is 
the dimer "internal energy" £' 2 ) = [U — 2j( 2 )], while 
the second term with the effective hopping rate J^ 2 ) = 
2J 2 /U gives rise to the dimer kinetic energy spanning 



the Bloch band e {2) {Q) = -2J^ C0S (Q) 0,0. Finally 



Hi, 



= V {2) rh 3 h ]±l 



W 
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H.c.) 



describes effective short-range interactions between the 
dimer and monomer, including a weak repulsive (or at- 
tractive, if U > 0) nearest-neighbor interaction = 
— 7J 2 /2U, and an exchange interaction with the rate 
W = 2J, the latter being responsible for the formation 
of off-site trimers as discussed below. 

In Fig.QJc) we plot the spectrum of the effective Hamil- 
tonian ([5]). It involves a two-body scattering contin- 
uum, with the energy given by the sum of the ener- 
gies of (asymptotically) free dimer and monomer, E c2 = 
[U-2JW] - 2 cos(Q) - 2 J cos(K-Q), and two bound 
states with energies E a \ and E a2 below and above E c2 . 
These dimer-monomer bound states are obtained using 
the Schrodinger equation for the two-body wavefunction 



^(Q, k) in momentum space, which leads to the integral 
equation 



*(Q,fc) = - 



I 

2^ 



dq 



Ul2 + Vcos(Q,q) + VsUQ,q) 
£( 2 ) +eW(q) + e(K - q) - E 



*(«,*), 
(6) 



where K = Q + k is the total quasi-momcntum, 
V cos (Q,q) = [2V^cos(q) - AJcos(K - g)] cos(Q) and 
V s in(Q,q) = V CO s(Q,q) with cos — > sin, while Ui2(— > oo) 
is an artificial dimer-monomer on-site interaction impos- 
ing the hard-core condition on Hamiltonian ([5]). Equa- 
tion ([6]) reduces to a non-linear equation for the energy 
E solving which (numerically) wc obtain E a i and E a2 
of Fig. [ljc). Comparison with Fig. [TJa') reveals good 
agreement with the exact solution for the three-body 
problem: the continuum spectra are practically identi- 
cal, while there are small but noticeable differences in 
the bound-state energies obtained from the exact and ef- 
fective models. These discrepancies are associated with 
the internal structure of the dimer [!, H, 0] , not accounted 
for by the effective model. We have verified that with in- 
creasing the interaction strength U, as the two particles 
forming the dimer become tighter co-localized, the dis- 
crepancies between the exact and effective models grad- 
ually disappear, as expected fH [Toj. 

For K = 0, ±7r, the bound states of Hamilto- 
nian ([5]) can be calculated exactly. To this end, we 
seek the eigenstates of H c s in coordinate basis, |\1>) = 
Z)ji^j a ^Uuh) \h,h), with the two-particle wavefunc- 
tion in the form $(ji,j 2 ) = e lK{jl+j2 ^ 2 e- lSK ^ (j) K {jr), 
where j r = ji — j 2l with ji and j 2 being the lattice 
positions of the monomer and dimer, and tan(<5ff) = 
tan (K/2)[J- jW]/[J+ J^]. For the relative coordinate 
wavefunction 0k (jr), imposing the hard-core condition 
4>k{0) — 0, we then obtain the difference relations 

Jk4>k{±2) + W k ^k(tI) + [E- V (2) ]<M±1) = 0, 
Jk[<PkUt + 1) + <Mir - 1)] + E<p K (j r ) = 0, (7) 



with \ Jr \ > 1, J K = y/ J 2 + J(2)2 + 2 JW> cos(K), 
W k = Wcos(K), and E = E - £^ 2 \ Equation © 
can be solved using the exponential ansatz <pK(jr > 

0) oc crjf" 1 and 4>ic{~3r) — = t0-R"(jV), which yields 

= ~Jk/[V^ =F Wk] for the symmetric (+) and an- 
tisymmetric (— ) wavefunction of the bound state {\chk\ < 

1) , with the corresponding energy E a in) = — + 

(a£°) a ]/a£>. 

It is now easy to see that without the exchange inter- 
action there would be no bound states (for any K) of 
the effective Hamiltonian ([5]). Indeed, this hypothetic 
(W — 0) problem is exactly solvable for all K, and 
for two hard-core bosons with nearest-neighbour inter- 
action there could be only one bound state when 
\o-k\ — \Jk/V^\ < 1 [1], which is not satisfied in the 



range of validity of the effective model [21(. In fact, 



Hamiltonian {5| can support two bound states because 



4 



exchange operators preclude the Bose-Fermi mapping 
[23 |. The effective nearest- neighbour interaction is, how- 
ever, responsible for the asymmetry in the binding ener- 
gies i?si(2) = -Eai(2)T2 Jk of the exchange-bound trimers 
below and above the continuum E C 2 [see Fig. QjU)]. With 
increasing the on-site interaction strength U, as — > 
0, we find that the binding energies approach asymptot- 
ically the universal limit Im^u^^ Eb\{2) = -F-J/2- In 
the general case of \K\ ^ 0, n, due to the exchange in- 
teraction, the eigenvalue problem H c s |W) = E\^) does 
not reduce to a closed expression for the relative coordi- 
nate wavefunction 4>K{jr), and the bound states cannot 
be calculated analytically. However, the universal limits 
for the binding energies Eb\(2) ~ ► -fJ/% remain valid for 
all K, since -> as |E7/J| -> 00. 

Discussion. We have found that, in a one- 
dimensional tight-binding periodic potential, three 
strongly interacting bosons can form two families of 
weakly-bound trimers with energies below and above 
the two-body scattering continuum of a single parti- 
cle (monomer) and an interaction-bound particle pair 
(dimer). Intuitively, these trimer states correspond to 
symmetric and antisymmetric states of a dimer and 
monomer at the neighboring lattice sites interacting with 
each other via an effective (particle) exchange interac- 
tion. In a two-site lattice, the splitting between the two 
trimer states would be equal to twice the exchange rate 
W = 2J. The infinite lattice, however, modifies the 
energy eigenstates: first a two-body continuum, corre- 
sponding to the dimer-monomer scattering states, opens 
up around the unperturbed energy of the dimer and 
monomer, [U — 2j( 2 '] ~ U, with the width given by the 
sum of the Bloch bands of the two, 4J + 4j( 2 ) > 4J 
(J*- 2 -* <C J for U 3> J). This continuum then pushes 
the modified symmetric and antisymmetric trimer states 
away from the edges of the continuum, resulting in their 
binding energies approaching I fJ/2 as U — > 00. 

The phenomena discussed above can be observed ex- 
perimentally with cold bosonic atoms in optical lattices 
|20|. The preparation of the trimers could be accom- 
plished by loading double-wells of a super-lattice with 
three atoms, followed by adiabatic conversion into a ho- 
mogeneous lattice. Using the optical Bragg spectroscopy 
(2^ | or lattice modulation rf spectroscopy followed by de- 
tection of the number of surviving bound pairs [H, one 
can then map out the bound trimer and unbound dimer 
energies of Fig. QJa) . 

Our effective model involving exchange interaction 
should be contrasted with the spin-^ fermionic case, for 
which the Bcthe ansatz solution 24j predicts no three- 
body bound state. Our results, therefore, apply to bosons 
only [25| . Studying larger number of bosons in a lattice 
might reveal other exotic bound states, while long-range 
interactions will certainly play an important role in the 
formation of bound states of three or more bosons in a 



lattice, as is the case in the two-body problem [§| . 
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